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We point out the importance of two-body final states of weak isosinglet neutral heavy leptons 
predicted in several models of new physics beyond the standard model. We concentrate on muon- 
type neutral heavy leptons L° with mass M < 2 GeV which can be searched for with increased 
sensitivity at a new round of neutrino experiments at CERN and Fermilab. Providing explicit decay 
rate formulae for the eev, e^iv, \x\iv, nfj,, pfj,, and ai[i final states, we use general scaling features to 
estimate sensitivity of L° searches in current and future experiments, emphasizing the importance 
of the -K/i decay mode. 



I. INTRODUCTION 

The LEP N u — 3 bound on the number of light, sequential-family neutrinos |lj tends to obscure the fact that 
relatively low-mass, neutral leptons, L , weakly mixed with one or more of e, \i, or r neutrinos, are still allowed. 
Grand unified theory models provide motivation, since some models may contain rather light neutral leptons |2j . Here 
we are specifically interested in L° with a mass M of order 1 GeV, and we assume for simplicity that a given massive 
neutrino mixes with only one light neutrino flavor. Experimental bounds exist ||, but it is timely to revisit them in 
light of the high intensity neutrino experiments coming on line or under development and construction. 



We especially emphasize the role of the L u 



z £^ (i^ = e^j /i^) decay mode in greatly increasing the sensitivity 



of searches for L° with M < 1 GeV, where phase-space strongly favors this two-body mode. For M > 1 GeV, the 
purely leptonic modes come into their own to provide clean signals for L° decays. 

II. HEAVY NEUTRAL LEPTON AND DECAYS 
A. Mixing 

In this study we are concerned with the search for heavy neutral leptons which are primarily isosinglets under weak 
SU(2)x,. These can mix with the light neutrinos, which are doublet members under SU(2)l. Whether the neutrinos are 
massless or massive but light, the mixing effect can be represented to a good approximation by replacing i>i h , i = e, fi, r, 
by 



5> 



aJ L 



(1) 



Here {v%) L and (£q) l are the left-handed components of the neutrinos and heavy neutral leptons, while (Mi) L are 
the combinations that appear in the charged and neutral weak currents (the weak interaction "eigenfields" ) . For 
simplicity we consider the case where each mixes with only one L°. If the neutrinos, z^, are massive, the mixing 
can be illustrated by a "seesaw-like" mass matrix B, 



M = 



where M. has eigenvalues 



A 

A M 



(2) 




(3) 



and corresponding eigenvectors 



1 




«-|^s_, ± ^=£_|. 



mq. + X 

When X/M « 1, L° = w+ w (A/A/, 1 - X 2 /2M 2 ) is associated with m+ « M, while 1/ = u_ (1 - A 2 /2Af, —X/M) 
belongs to m„ = to_ w X 2 /M. If, as in some specific models, A were related to charged lepton or quark masses, the 
mixing parameter U ~ X/M could be determined by the heavy neutral lepton mass alone. 

In the work we present here, U and M are treated as independent parameters. Since our treatment is perturbative, 
the mixing indicated in Fig. |l| can be characterized by a mixing mass parameter A, while the L° and v lines represent 
propagators or a propagator and an external wave function. For example, treating L° decay, one has factors 

a4 ± ^ [ W) = T7 —Ulo(p) (5) 

where C/x,o(p) is the L° momentum space wave function and m v /M « 1. The parameters A and M are to be 
determined, or bounded, by experiment. 

FIG. 1. L° mixing. 

Figures ^ and ^| show the Feynman diagrams for the charged and neutral current contributions to leptonic LP a decay 
in lowest order in £/j 's and weak couplings.^ For a final state where li ^ £j and both leptons are observed (which is 
required in the analysis presented here), only the charged current processes, Fig. ^| contribute. Assuming the charges 
of the final leptons are determined, say £~ and £f, then there is a single diagram with factor Ui a . If the flavors, but 
not the charges, of the leptons are determined, one has a second mode and another diagram with £~ at the upper 
vertex and if at the lower one. The decay involves two terms, one with a factor \Ui a \ 2 and one with |J7j a | 2 - For the 
en case, T e ^ + = r^~ e+ , so r e ^ = {\U e \ 2 + \U^\ 2 )T^ e+ . 

FIG. 2. L° charged current decay. 



FIG. 3. L neutral current decay. 

In the case that li = £j, the neutral current diagrams are involved, and each Vi can contribute for a fixed j and fixed 
a. There is only one charged current contribution, proportional to U a j, in this instance. The diagram where i = j 
in Fig. [|, also proportional to U a j, interferes with that of Fig. ^, while the other two terms with i ^ j contribute 
incoherently to the partial rate for L° a — ► £j£j+ missing neutrinos. We give the expressions for the muon-type L° in 
Appendix |A|. The numerical results and parameter bounds that we present will be for the situation where only muon 
neutrinos mix with a given heavy neutral lepton. The results can be generalized straightforwardly to the mixed cases. 

While on the topic of mixing, it is worth offering the reminder that heavy neutral leptons with masses below 2 GeV 
will live long enough to escape the LEP and SLC detectors before decaying and thus they are disguised within the 
N v = 3 light-neutrino number measurements [Q. The right handed components of L° do not couple to Z°, while the 
left handed components couple through the mixing displayed in Eq. |l| and in Fig. [I], which ensures that the "missing" 

decay partial width, which goes like I | 2 j gives the usual result as long as L°'s are light enough to be 

i=l 

treated to the accuracy of present measurements as massless particles in the final state. 



x For definiteness we diagram the decay as if the source beam were Dirac NHL's. If the beam is sign-selected (like that of the 
NuTeV experiment, for example), then this case applies directly. If the neutral lepton source is from a beam dump (as in the 
DONUT experiment, for example), then there is a democratic mixture of neutrino and antineutrino; and the conjugate L°'s 
and consequent conjugate lepton final states are all present. For Majorana-type NHL, each charge state can occur in any decay 
regardless of the nature of the beam preparation, and the diagram with the lepton switched in Fig. Ill also occurs. 
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B. Decays 



In this section, we present the partial decay widths of L° to channels of relevance for our study. For defhritencss we 
restrict our attention to a "ix-type" L° that mixes only with the weak muon neutrino with mixing strength denoted 
by U. 



1. Leptonic Decays 

Referring to Fig. 2, we see that the relevant decays are 

L° /i-e+^ (6) 



and 

P P "P 



L° - M + M -^ . (7) 



For convenience, we divide out a common factor 



g\m\ t2 



192tt 

which happens to be the all-neutrino decay rate summed over flavor 



K = TT^ZT^ . ( 8 ) 



K= Y, ? (1% ^ u^m) , (9) 



e,/i,r 

and to the usual //-decay width formula when the electron mass is neglected, up to the factor U 2 . Differential and 
integrated partial width formulas are summarized in Appendix |a|. Integrated partial widths, r(X" — * [i + e~ve)/K 
and T(L° — > jtx + pTv^)l K are shown as a function of M, the L° mass, in Fig. [|. Clearly the three-body leptonic 
modes come into their own above M = 1 GeV, but the hadronic modes summarized next also play a major role in 
detection of L° decay. 

FIG. 4. The important leptonic and hadronic decay widths, T/K, as a function of M in GeV. The solid curve is the efiv 
mode, dotted is fifiv, dashed is 7T/x, dot-dashed is pfj,, and dot-dot-dashed is ai/x. Note that the 7r/x peak is more than 7 times 
higher than the e\w peak. 



2. Exclusive hadronic decays 

The exclusive decays 

Ll->^-H+, (10) 
where H = ir,p or oi, are large modes that can be reliably calculated within our framework. As we will explain in 



Sec. IIIB, only partial widths are needed for experimental searches since U is already constrained to be small. The 
signatures are clean, and the L° detection sensitivity increases dramatically when these search modes are included, as 
we discuss in the following section. The isospin-related neutral modes, L° — > v^H°, are not as useful experimentally, 
and are not explicitly included in the discussion. 

The L° — > /i~7T + partial width is fixed by 7r + — > [i^v^ while T \L° — > n~ p + ) and T — * /x^a^J are fixed by 
t + — > and r + — * ^ r decay. The ratios r(L° — * zx~ H + )/K are shown in Fig. |as a function of M. The special 
prominence of the L° — ► pTn^ mode for low mass L° is evident. Regarding the interpretation of the figure, note that 
two-body decays appear to rise, peak and fall as a function of M because of our normalization to the 3-neutrino decay 
width K. Two body decays are proportional to G 2 F M Z f^, where fu = / w in the pion case, for example. The peak 
in the curves for the 7r + ,p + and final state partial rates is then an artifact of the M~ 2 from our normalization 
multiplied by the rapidly rising two-body phase space. 

As an aside, we note that two-body decays allow the potential discrimination between electron-type and muon-type 
NHL through the observation of L° — * fj,~H + vs L\ — > e~H + . 

Next we turn to a consideration of the issue of sensitivity to detection in experimental searches, focussing on 
those aspects that follow from basic scaling considerations and are independent of details of specific experiments, 
which must be handled by thorough Monte Carlo analysis. 
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III. ESTIMATES OF SENSITIVITY 



A. Production of L° from Meson Decay 

Light neutral heavy leptons could be produced in decays of any charged stable meson; K ± , D^, and provide the 
most experimental sensitivity because they can be produced in enormous quantity in the pN collisions that generate 
neutrino beams and because CKM suppression of the two-body decay is modest(i4T ± , Z? ± ) or absent (Dg). The 
two-body decay rate can be related directly to the [iv decay rate 



T(H -» M L°) = U 2 T(H -> ^) VP (M/M H ,M^/M), 



(11) 



where H = (K ± , D ± , Dg) , U 2 is the mixing factor, Mjj is the meson mass, M is the L° mass, is the muon mass, 
and 



(l + y 2 )~ X 2 (l-y 2 Y 



l~x 2 (l-y) 2 l-x 2 (l + y) 2 



y 2 (1 — x 2 y 2 Y 



(12) 



is a kinematic factor. Note that for M 2 /M 2 > 1, ■q P (M/M H ,M ll /M) -> (M/Af M ) 2 (l - (M/M H ) 2 ) 2 , demonstrating 
that two-body decays to L° can be considerably enhanced by the lifting of helicity suppression. 



B. Experimental Detection of L° 

The simplest scheme for a detector is a low mass decay space immediately in front of a neutrino detector 
instrumented with tracking chambers to permit reconstruction of a possible decay vertex. By removing as much 
mass as possible, for example through the use of helium-filled bags, backgrounds from conventional v^N interactions 
can be minimized. The neutrino detector downstream of the low-mass region can be used to identify muons, pions, 
and electrons and measure their energies, thus providing sensitivity to the largest decay modes of L° with M < 2 
GeV, L® — ► e/i, fifi, ivfi. While a dedicated instrumented low-mass detector is optimal, it is not necessary, provided 
the density of material in front of the neutrino detector is not too high. 

If a decay space is constructed to have a length A along the direction of the beam and a width much wider than 
the beam, the probability of a being observed in the detector can be expressed as 



p2 i M,U* ) =r A ' 1 ^e-->*^%$, (13) 

Jz-A/2 7PCTL0 1 L o 

where Tjo',T t fo are the detected and total decay widths, z is the distance of the decay position from the LP production 
point, Z is the distance from the production point to the center of the decay channel, (3 is the L° speed in units with 

c = 1, 7 = l/^/l — f3 2 , ^Pct^o is the L° mean decay length, and erf is the average L° detection efficiency. Mixing 
angles of interest are sufficiently small and detectors are configured such that A -C Z 7/3ct £ o . In this case, Eq. [l3| 
simplifies to 

AT^dct 
1 T° r° 

Pn»{M,U*) = ——-£e I f (14) 

1PCT L 1 L0 



pdct 

A._^/° 

7/3 he D 



(15) 



Since r^o Ty * = h, the observation probability calculation requires only knowledge of the partial decay widths of the 
channels being searched for, rfo*. As these partial widths can be reliably calculated, there should be little theoretical 
uncertainty in estimates of search sensitivity. 
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C. Sensitivity in (U 2 ,M) Plane 



The number of produced L° will be proportional to the total number of protons on target, iVpoT, and a factor 
U 2 M 2 from Eq. [0], assuming M 2 -C M 2 -C M#. Using the simple efj, decay mode as an example, the detected 
partial width will be proportional to U 2 M 5 . The number of decays in the decay region is proportional to this partial 
width, to the length of the decay space A, and to an extra factor of M/E from time dilation, with E the L° energy. 
Combining all effects together, the number of observed L° has the dependence 

N obs ^ %A[/4M 8 (16) 
M E 

If an experiment performs a search and observes a statistically significant L° signal, then the experiment can determine 
both M, from the two-body semihadronic decays, and U 2 from a more detailed development of Eq. [l6| (see Appendix 
|B|). On the other hand, if no candidates are observed, then the experiment can exclude a region in the [U 2 ,M) plane 



that from Eq. 16 will be bounded by a curve of the form U M = constant. The minimum mixing factor sensitivity 
for fixed M will be proportional to y i?/(7VpoT A) for an experiment that suffers no background. Greater sensitivity 
follows from increasing iVpoT and A or decreasing E, but the gain is slow because of the square-root factor. Adding 
the hadronic decay channels presented here is another way to improve the search limits. 



D. Estimates of Sensitivity for Current and Future Neutrino Experiments 

This section summarizes predictions for two rather different experiments: NuTeV p"2] ], a high energy deep inelastic 
scattering experiment with its neutrino detector located far from the production target, and DONUT Jl3| , a high 
energy beam dump experiment with its detector 35 m from the neutrino production target. We also comment on 
prospects in lower energy experiments. A more detailed prescription for our estimates is given in Appendix 0. 



1. The NuTeV Experiment 

NuTeV has installed an instrumented decay channel to search for NHL [^) and is currently taking data. The 40 m 
long NHL decay channel is located approximately 1200 m from the decay position for charmed mesons and 900-1200 
m for charged kaons. NuTeV may receive an integrated intensity of up to 6 x 10 18 protons-on-target. Neutrino 
interactions occur at a rate of approximately 20/10 13 POT. The ratio of kaons to pions in the NuTeV beam is about 
0.4 resulting in v e /v^ interaction ratio of about 2.3% in the detector. Contributions from charmed meson decay 
increase the v e rate by approximately 1% of itself. The Lab E neutrino detector used by NuTeV has a fiducial length 
of approximately 15 m and a mean density of 4.2 g/cm 3 . 

Figure || shows the estimated sensitivity plot in JJ\ vs M for NHL produced from decay in NuTeV. The main 
model sensitivity comes from the assumption for the mean K + energy, which enters both in the lifetime calculation 
of the L° and the interaction probability for the neutrinos used in the normalization sample. 

FIG. 5. Estimate of sensitivity as a function of mass M (in GeV) to NHL produced from kaon decay in NuTeV. The vertical 
axis represents the minimum mixing parameter \U\ 2 the experiment would be sensitive to from kaon decays alone assuming no 
candidate events were observed (and no background events were expected) in an exposure to 6 x 10 18 protons-on-target. The 
dashed, solid, and dot-dashed curves are for kaon energies of 100, 150, and 200 GeV, respectively, and assume the search is 
performed only using the L° — > e/i decay mode. The dotted curve illustrates the sensitivity that could be gained for 150 GeV 
kaon energy if the 7r/x, p/i and ain modes were added to the e\i mode in the search. 

Figure || shows a similar plot for NHL produced from Dg decay. The main model sensitivity in this result comes 
from the assumption for the mean D§ energy. 

FIG. 6. Estimate of sensitivity as a function of mass M (in GeV) to NHL produced from Ds decay in NuTeV. The vertical 
axis represents the minimum mixing parameter \U\ 2 the experiment would be sensitive to from kaon decays alone assuming the 
no candidate events were observed (and no background events were expected) in an exposure to 6 x 10 18 protons-on-target. 
The dashed, solid, and dot-dashed curves are for charmed hadron energies of 50, 100, and 200 GeV, respectively, and assume 
the search is performed only using the — > efi decay mode. The dotted curve illustrates the sensitivity that could be gained 
for 100 GeV Ds energy if the tt/j,, p/i and oi/x modes were added to the e\i mode in the search. 
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Figure [7| shows the estimated sensitivity for all modes combined using either the e/x decay channel by itself or the 
e/Lt and wfx channels combined. Note that the U A dependence of the detection rate for NHL implies that uncertainties 
in acceptances, branching fractions, etc., which only enter as ratios, affect results only as square roots in determining 
U 2 . The dependence of the result on statistics, assuming no background, is similarly proportional to 1/V-^Vpot- 

FIG. 7. Estimate of combined sensitivity to NHL from K, D, and Ds decays in NuTeV as a function of NHL mass M (in 
GeV). The vertical axis represents the minimum mixing parameter \U\ 2 the experiment would be sensitive to from all meson 
decays combined assuming no candidate events were observed (and no background events were expected). The curves assume 
6 x 10 18 POT, a mean kaon energy of 150 GeV, and a mean charmed meson energy of 100 GeV. The solid curve shows the 
expected result assuming the search were performed only using the L° — » e/x decay mode. The dotted curve illustrates the 
sensitivity that could be gained by adding the tt/j,, p/i, and a±fi modes to the search. 

2. The DONUT Experiment 

DONUT is a hybrid emulsion spectrometer detector sited approximately 35 m from a beam dump target in Fer- 
milab's 800 GeV proton beam. Active and passive shielding eliminate essentially all neutrinos produced by pion and 
kaon decay, leaving a mixed beam of v e , v^, and v r from charmed hadron decay. The experiment's primary goal is to 
detect charged current interactions of v T with nucleons in the emulsion target. DONUT may receive an exposure of 
up to 2 x 10 18 protons-on-target. 

DONUT's proximity to the production target greatly enhances the flux of neutrinos, and possible NHL, produced 
from D ± and D s decay. Only a few thousand charged current interactions in the emulsion detector are expected, 
but these will essentially all originate from charm decay. Figure || shows an estimate to sensitivity |l4| for NHL in 
DONUT, assuming the experiment could instrument a 5 m decay space in front of their emulsion target, and that 
the experiment receives an exposure of 2 x 10 18 POT. The estimate is comparable to that for NuTeV, and possibly 
better at high L° mass. 

FIG. 8. Estimate of combined sensitivity to NHL from D and Ds decays in DONUT as a function of NHL mass M (in 
GeV). The vertical axis represents the minimum mixing parameter \U\ 2 the experiment would be sensitive to from all meson 
decays combined assuming the no candidate events were observed (and no background events were expected). The curves 
assume 2 x 10 18 POT and a mean charmed meson energy of 100 GeV. The solid curve shows the expected result assuming the 
search were performed only using the L° — > e/i decay mode. The dotted curve illustrates the sensitivity that could be gained 
by adding the nfi, pp, and a\\i modes to the search. 

3. Other Experiments 

We have used the NuTeV and DONUT experiments as specific examples in calculating sensitivity to NHL; however, 
our formulas can easily be applied in other cases. CHORUS fl6|| , also a hybrid emulsion-spectrometer, and NOMAD 
p5[ , a low mass high resolution spectrometer, are running in a low energy horn beam at CERN to search for v T 
produced from — > v T oscillations. Two additional — > v T oscillation experiments, COSMOS jl7| and MINOS 
p8| at Fermilab, will begin taking data around 2001. All of these experiments can search for NHL produced from 
(primarily) kaon decay because of their lower energy beams. We note that searches in this mass regime benefit 
appreciably from an ability to detect the ir/j decay mode. 

IV. SUMMARY 

As just discussed, a number of Fermilab and CERN neutrino oscillation experiments are currently running or will 
soon come on-line, and we have shown that a simple and direct method to expand the search for light (M < 2 GeV) 
neutral heavy leptons (NHL) is suitable for all of these experiments. The present lower limits on lifetimes for NHL's 
in this mass range means that all experiments satisfy the criterion that A << z << jPct^o where A is the fiducial 
length, z is the source-to-detector distance and 7/3cr is the decay length. Therefore the simple criterion (|l4|) applies. 
Only theoretical values of the partial widths into search modes are relevant, and these can be reliably calculated in 
terms of the mass and mixing parameters. 

In addition to common features that lend themselves to a clean analysis, we have shown how important the two- 
body decay modes are, especially the irpi mode, in achieving improvements in sensitivity. We found significant gains 



G 



over most of the mass range M < 2 GeV, which means discovery reach extends to smaller mixing values than currently 
reported in the literature. 

Individual experiments will use sophisticated Monte Carlo simulations to properly account for details such as 
efficiencies and fluxes; but the simple and general formalism we have presented applies to all experiments. We believe 
our results provide a broad, useful framework for expanding searches for neutral heavy leptons. 
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APPENDIX A: EXPLICIT DECAY RATE FORMULAE 



In this appendix we collect formulas for partial decay widths to which we refer in the text. For completeness we 
also include the differential forms before the final integrations, including effects of L° polarization in the e/i decay 
channel. All formulas refer to the L° rest frame. 



1. Leptonic Rates 

Referring to Fig. || and identifying jjT = ii and e + = £j, with x^ = 2E IXr /M, x e = 2E e /M, and x m = M^/M, we 
have 



d 2 r ( M e+) q2 M 5 



dx p dx I 



16tt 3 



|E^| 2 M1 -X e - X 2 m )} , 



(Al) 



where M is the L° mass, M M is the muon mass and E e and are the positron and muon energies. We have neglected 
the electron mass in writing Eq. Al. Saving the x^ integration until last, and setting m e = in the phase-space 
treatment, the single differential decay form is given by the expression 



dT^ e ) 

dxu, 



l07T J 



2 i 



— X,, 

3 



(A2) 



Integrating over Xp, from 2x m < x^ < 1 + x^ , we get the result familiar from muon decay, 
For completeness we include the corresponding expressions for the L° e —>■ e~ n + Vp final state. 



(A3) 



where again the electron mass is neglected. Integrating over the electron energy first, we obtain 



dx,, 



G%M 5 
16tt 3 



\U e \ 2 x^(l- x^ + x 2 n ) 



2 _ A T '2 



(A4) 



(A5) 



Comparing Eq.'s A2 and A4, we see the interesting feature that, in principle, the difference in the two distributions 
could distinguish — > \iTe^v e from E\ — ► e~ [i^v^ even if the charge states were not determined directly. This is 

true even if x m << 1, though distinguishing between distrib utio ns e+ ) 
be difficult. Integrating over the muon energy, one finds Eq. A3. 



?, - la;?, and r( e ~ x 2 ^ — x^ may 



3j" 
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If the L° is produced through meson decay, it will likely be highly polarized. Correlations between the polarization 
vector and the e and [i momentum vectors could produce substantial effects in detection efficiency for the L° — > [i~e + v e 
final state. To permit study of these effects, we give the following triply differential decay formulas: 



d 3 T ( P «,+) GlM 5 lTT l2 , 2 



dx^dx e dcos8, 



32tt 3 



(A6) 



where 9 e is the angle between the final e + direction and the polarization direction of the decaying L°,and m e is 
neglected as before; or, alternatively, 



|[/ M | 2 iE e (l - x 2 m - x e ){\ - cos0 e cos0„), 



dx^dxedcosO^ 



32tt 3 



where 



cos0» — 



(2 x^ x e ^j x „ x 



4xi 



-v / %ft ^^rn 



(A7) 



(A8) 



and Oft is the angle between the final /i - direction and the polarization direction of the decaying L . 

Turning to the case where there are two muons in the final state and the charged and neutral current contributions 
interfere, we write the double differential decay formula as 



d 2 T^ ^ 
dx+dx- 



C 2 M 5 

-^-^|ig 2 [x + (l - x+)(a + bf + x_(l - x_)(a - bf 

„2 to „ ™ \/"„2 ,2s 



(A9) 



+2<(2 -x--x+){a z -b z )] , 
where x± — 2E„± /M, a — (3/2—2 sin 2 6y/) an d 6 = 1/2, and sin 2 9yy = 0.224. Integrating over x + , one obtains 



dT^ * ) G 2 Af 



.^{Ml -*_)(<, -ft) 5 



16^ 3 1 Ml 
+2a&(2 - x_)(a 2 - 6 2 )] {x u _ - xL) 

+ha + &) 2 - 2, 2 > 2 - 6 2 )] [(xir-(xtr 



(A10) 



The expressions for x"l and x_ in Eq. A10 read: 



(l+x^-x-) 



(AH) 



One obtains the corresponding expression for dT ^ x - — by the replacement Integrating the expression |A10| 

over X-, one obtains the partial width formula for L° — > namely 

r(M - M +) = | Cl [ (: _ 14 ^ _ ^ _ 12^)^1— 4^ 



1927T 3 



+12^(^1 - + 4C 2 ^(2 + 10a; 2 , - Ux^)y/l - 4a* 
+6x^(1-2^ + 240^} , 



where 



L = log 



l-3a; 2 „-(l-x 2 ra ) v /l-4a : 2 

^(1 + ^1- 4 ^m) 



(A12) 



(A13) 



The coefficients Ci and C2 in Eq. A12 are 
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Ci = - [1 + 4 sin 2 d w + 8 sin 4 W ] 
C 2 = - sin 2 9 W [1 + 2 sin 2 W ] • 



For completeness we give the corresponding L° — * e + e ^ partial width 

r (e-e+) = G|MJ 2 1 r 1 _ Asin 2 



'W 



sin 



(A14) 
(A15) 

(A16) 



2. Two— body Semihadronic Decay Formulas 

The three decay modes of interest in this work are L° — > tt + \iT , p + p~ and a^/x _ , which involve the decay constants 
f % , <?p, and g . We determine g p and <? a from the partial widths for r — * pv T and t — * a\v T using the lowest order 
diagrams. As justification, we note that the t — ► itv T partial width calculated with fusing the measured value of / w 
agrees with the experiment to within a few percent. The decay formulas in the rest frame of L° read, in the narrow 
p, a\ width approximation 



dQ 



ff-^Gl\U^M\/S(M, M7T ,M,) 



K 
M 2 



Ml 

M 2 



Ml 
M 2 



(A17) 



dn 



1 + 



M 2 \ Mi 



Mt 



M 2 / M 2 M 4 



M 2 



(A18) 



where 



S{M,M H ,M tl ) 



I ( Mh My. 

M M 



1 - 



M H M p 
M M 



(A19) 



and dT^ aie ^/dil is obtained from Eq. A18 with th e rep lace ment s g p — » g a and M p — » M a . The integrated partial 
widths are obviously obtained by multiplying Eq.'s Al7| and A18 by \-k._ 



The parameters g p and g a are determined from the r partial widths |8J to be 

r(r -> pv T ) 
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and similarly 



fc 0S ^(l-^)( 
(0.102 GeV 2 ) 2 , 

5 2 = (0.128 GeV 2 ) 2 . 



1 



ITT 



(A20) 



(A21) 



The branching fraction for t — ► ir n + ir + > neutrals +i/ T was used for the a\ fraction. A recent chiral dynamics 
analysis of mesons M yields 



and 



g p = 0.104 GeV 2 



g a = 0.136 GeV 2 



(A22) 



(A23) 



in reasonable agreement with the values in Eq.'s A20 and A21 







APPENDIX B: DETAILED ESTIMATES OF SENSITIVITY 



A useful way to estimate sensitivity to NHL couplings is to normalize to neutrino interactions in the detector 
produced by the same parent meson as that which produces the NHL. This technique relies less on absolute calculations 
from a beam monte carlo and allows one to identify important model dependencies. 

This section will develop estimates of sensitivity to NHL produced from K^, D^, and D$ two-body meson decays. 
At the end we comment on NHL search possibilities in all current and future accelerator based neutrino experiments 
that we are aware of. 



1. ^i-type NHL from Kaons 



We define the observable 



„ Detected L®, in Channel from K + Decays 

R (M U ) = - (Bl) 

Detected in Detector from K + Decays 

- R P (M, U 2 )R B (M, U 2 )R d (M, U 2 ), 

with Rp(M, U 2 ), Rb(M, U 2 ), and Rd(M, U 2 ) defined as the L° production ratio, beam transport ratio, and detection 
ratio, respectively. The number of detected L° will be a strong function of L° mass M and mixing factor U 2 . The 
number of ordinary detected from K + decays can be inferred from the energy spectrum of neutrino interactions 
in the neutrino detector. By formulating the search as a measurement of the ratio Rk{M, U 2 ) 1 one lessens sensitivity 
to absolute normalization of the neutrino beam and detector acceptance. To establish the limits of experimental 
sensitivity, we will consider the null case, where no L° candidates are observed for a given exposure in a neutrino 
beam. In this case, the upper 90% confidence level limit sensitivity for Rk(M,U 2 ) follows, assuming no observed 
events and no background, from Poisson statistics, 

Rk(M, U 2 ) < ^ , (B2) 

with tk-k the fraction of events from kaons, n v the number of interactions in the detector per incident proton, 
and -/VpQrp the total number of protons on the production target. This section will derive the functional dependence 
of Rk(M, U 2 ) on M and U 2 that will allow limits to be placed on the (M, U 2 ) plane from a null result search. 



a. Production Factor Rp(M,U 2 ) 



Rp(M, U 2 ) is the ratio of produced NHL to produced from kaons M: 

U 2 M 2 (i " , M M„ 



Mi 



where, as before, U is the — L° mixing strength, rjp is defined in Eq. [l2| and M, M^, and Mk are the NHL, muon, 
and charged kaon masses, respectively. If the muon mass is neglected, the phase space function r]p(x,y) is unity. 



b. Beam Factor R B (M, U 2 ) 

Rb(M,U 2 ) is the relative acceptance for NHL vs neutrinos due to the beamline. An experiment will generally 
run a detailed Monte Carlo simulation to obtain this factor, but at high energies, E L o M, one would expect only 
modest acceptance differences. Accordingly, we assume for estimation purposes that 

L l 

R B (M,U 2 ) = ^-^ 1.0. (B4) 
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c. Detection Factor R D (M,U 2 ) 



Rd(M,U 2 ) is the relative detection probability for NHL vs v^. This can be expressed as a ratio of detection 
probabilities, 



R D (M,U 2 )= puK ■ (B5) 



L° 

For long NHL lifetime, the NHL detection probability P D "(M, U 2 ) can be written 

1 TO r° 

F D "(M,U') = — -±e L D \ (B6) 

where Fro t ,r t 7 ° t arc the detected and total decay widths, A is the length of the decay channel, ^ficr^o is the L° decay 

L° 

length, and Ejf is the mean L° detection efficiency, assumed to be 1.0. Choosing to observe only the e/x mode for 

L° 

simplicity, P D ''(M, U 2 ) can be re-expressed as 

r° TI 2 M e A r° 
i#*(M, U 2 ) = , s D ^ D (M,/M), (B7) 

with r]D ("m 1 ) a threshold factor in the L° decay that can be read from Eq. A3 for the e/i decay. 

The neutrino interaction probability is simply 

P» D « (M, U 2 ) = N A (pt) a' u (E v ) e v « , (B8) 

where a' v (E v ) — 0.68 x 10 - 38 (E v ) cm 2 is the neutrino-nucleon cross section, Na (pt) is the target thickness in units 
of nucleons/cm, and e^f is the target detection efficiency, which is usually close to 1.0. 



d. Sensitivity Formula 



The final prediction for Rk(M, U 2 ) is of the form 



(B9) 



with C containing all experiment-specific information, 



C = a- 



and a a constant, 



E L0 ){E v )e v «{pt) 



K f 1 - 5) ° r * N *< 



(BIO) 



= 2.7 x 10 



16 cm- 3 E 



^GeV 



For fixed M, it follows from combining equation B9 with equation B2 that 



Uk < M ~ 




M 2 
Ml 



2.3 



(Bll) 

This formula demonstrates the qualitative feature of any NHL limit based on a search for decay vertices of NHL 
produced from meson decay. Sensitivity to the mixing U\ is proportional to M~ 4 up to masses near the kinematic 
limit. The M -4 behavior originates from the U 4 M S dependence of Rk(M,U 2 ), with U 2 M 2 arising from NHL 
production, U 2 M 5 from NHL decay, and an extra factor of M from time dilation. Because of the U 4 dependence of 
Rk{M, U 2 ), limits on U 2 only improve as the square root of the total integrated protons on target. The square root 
also insures that other experimental effects, which enter as ratios in C, only weakly affect the sensitivity estimate. 
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2. /i-type NHL from D+ Mesons 

The calculation proceeds in a manner similar to the kaon case, except this time a more convenient normalization is 
relative to electron-neutrinos from charmed mesons: 

„ Detected lP u in Channel from D + Decays 

Rr(M U 1 (B12) 

' Detected v e in Detector from D + Decays 

= R' P (M, U 2 )R' B (M, U 2 )R' D (M, U 2 ). 

This choice is motivated by the fact that charm decays produce a negligible fraction of the total v^, but on the order 
of 1% of the v e in a high energy beam. The 90% C.L. Poisson sensitivity to Rc(M, U 2 ) can be written as 

Rc(M, U 2 ) < , (B13) 

where f„ e is the fraction of neutrino events which are v e and ] Ve D is the fraction of interacting electron neutrinos 
from D + decay. Note that f Ve D — > 1 in an ideal beam-dump experiment, whereas f Vc D — 0.01 in a conventional high 
energy experiment. 

For the production rate, one can assume that all conventional v e are produced from three-body scmilcptonic decays 
whereas the NHL production is dominantly two-body, 

RAM, U ) - [I'-jr) B{D+ ^^ x y 

In this formula the muon mass is neglected compared to the D + mass Mjj, and the ratio B(D + — > n + v^) / B(D + — » 
H + X) corrects for the use of three body decays of the D + as the dominant source of v e . The D + muonic rate can 

be calculated assuming a D meson decay constant fjj |l(| via B(D + — > n + v^) = (3.2 x 10~ 4 ) ( 200 M e y ) ' anc ^ ^ e 
semi-muonic rate B(D + — ► n + X) is measured as B(D + — > n + X) = 0.172. 
For the beam transport factor, we assume again a value of 

R' B (M,U 2 ) = 1.0. (B15) 

This number may be less than 1.0 owing to the higher px given to the L° relative to the three-body decay neutrino. 

The detection factor is modified by the softer D + spectrum, relative to kaons, and the somewhat higher L° energy 
expected relative to the v e from three-body decay. Again, take 



P D "(M, U 2 ) 



r' d (m,u 2 ) = ^ d vj;r ; > ( Bie ) 

D 



with 



and 

P'» K (M, U 2 ) = (pt) u' v (E„ cD ) N A e^ D (B18) 
The final expression for Rc(M, U 2 ) is 

R c (M,U 2 ) = C , U i M a i [ l~-^-j , (B19) 

with 

'A 



C 
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and 



M , 




, B{D+ -+ y+X) 
B{D+ 



4.6 x 10 13 cm _3 g -1 GeV 



-6 



For fixed M, it follows that 



[/£ < M -4 




M 2 



2.3 



(B20) 



Sensitivity extends to much lower mixing angle than the kaon case, essentially because the M 5 behavior of the 
decay rate overwhelms the ^1 — jtfirj phase space factor, Charmed particles are also typically produced with lower 
energy in a hadron collision than kaons. This implies that L° from charmed hadron decays will have smaller time 
dilation factors, E' L „ /M, and so the probability that an L° will decay in front of the detector is larger if the L° 
originates from a D + decay than it is for L° from K + decay. 



3. M-type NHL from DJ Mesons 

The limit on U 2 from this source follows directly from the result for D + mesons since one can immediately relate 
the number of NHL from Dg decay to the number from D + decay via 



B(D+ - n+LO) ~ V 2 d M D {\ 



M 2 \2 f2 



(B21) 



and the ratio of produced cs mesons to cd mesons pT| , 

N(Dp 



r D s D 



0.5. 



N(D+) 

It follows that a limit on NHL produced from Dg is directly related to that from D + by 



(B22) 



\ 



V 2 M D (1 

cd 



M 2 \2 



f. 



1 



V c 2 s M Ds {l 



M 2 



) 2 fL r DsD 



(B23) 



Dg contribute more sensitivity to the rate than D + owing to the Cabibbo favored annihilation diagram in the former's 
decay. This estimate's main uncertainties are the ratio of Dg to D + production, rjj s £>, which is taken to be 1/2, and 
the ratio of pseudoscalar decay constants, /d//d s , taken to be 1.0. 
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